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ABSTRACT 


The transverse forces acting on smooth and sand-roughened circular 
cylinders immersed in an harmonically oscillating flow have been 
measured using a recently constructed U-shaped water tunnel. 

The lift coefficients were found to depend on Keulegan-Carpenter 
number but independent of the Reynolds number within the range of 
relative roughness and Reynolds numbers tested. The results have shown 
that the lift coefficients for smooth cylinders at low Reynolds numbers 
are nearly identical with those obtained for rough cylinders at very 
high Reynolds numbers at the corresponding amplitude ratios. The results 
have also shown that the transverse force is a significant fraction of 
the in-line force and must be taken into consideration in the design of 


structures. 
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NOMENCLATURE 
Virtual amplitude of the motion 
Amplitude of the motion 
Maximum acceleration 
Maximum transverse force coefficient, also used as C, (max) 
Diameter of the test cylinder 
Instantaneous total force acting on the test cylinder 
Frequency ratio, Vg WL 
Frequency of the first harmonic of the transverse force 
Measured force 
Gravitational acceleration 
Distance between pressure taps and mean water level 
Keulegan-Carpenter number 
Length of the test cylinder 
Reynolds number 
Sand roughness particle size 
Distance between the pressure taps 
Period of oscillations 
Time 
Maximum velocity 
Instantaneous velocity 
Width of the test section 
Frequency parameter 
Specific weight of fluid 
Differential pressure 
Phase angle 
Fluid kinematic viscosity 


Fluid density 
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I. INTRODUCTION 


It is a well known fact that flow about bluff bodies gives rise to 
separation and vortex shedding. Aside from the fact that the formation, 
growth, and shedding of vortices produce a drag force on the body in 
the direction of mean flow (hereafter referred to as the in-line force), 
the alternate shedding of the vortices also results in a time-dependent 
transverse force. It is of both fluid-dynamical and practical importance 
to evaluate the amplitude and the frequency of all the harmonics of the 
transverse force. The primary reasons for this are that the transverse 
force may set the body in motion resulting in a so-called hydroelastic 
oscillation and fatigue, and that the vectorial sum of the in-line and 
transverse force may be the most important force for the design of a 
structure. 

Evidently, the transverse force depends on the shedding characteris- 
tics of vortices. These in turn depend on the characteristics of the 
ambient flow (e.g. steady flow, harmonic flow, wavy flow, etc.), on the 
Shape and orientation of the body relative to the flow, on the surface 
characteristics of the body (smooth or rough), and on the parameters 
characterizing the fluid-body interaction (Reynolds number, Keulegan- 
Carpenter number, wave height to depth ratio, etc.). 

Of the scores of papers written on the fluid loading on structures, 
none seems to have systematically treated the effect of surface roughness 
on transverse force, particularly at high Reynolds numbers. Previous 
Studies, such as those carried out by Chang [1], Bidde [2], Wiegel and 


Delmonte [3], Mercier[4], Sarpkaya and Tuter [5], Jones [6], Isaacson [7], 
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and Sarpkaya [8] dealt with smooth cylinders at relatively low Reynolds 
numbers. Only recently, Collins [9] has conducted a series of experi- 
ments with sand-roughened cylinders at fairly high Reynolds numbers and 
demonstrated the importance of surface conditions. The data obtained 

by Collins have been extended in the present work to cover a wider range 
of relative roughnesses and Reynolds numbers. Thus, the primary pur- 
poses of this continuing investigation were to furnish data, obtained 
under carefully controlled laboratory conditions, about the transverse 
force acting on sand-roughened cylinders with relative roughnesses from 
1/800 to 1/50 in harmonically oscillating fluid at Reynolds numbers up 
to about 1,000,000 and to attempt to identify the physical mechanisms 
responsible for the correlation or scatter of the transverse force co- 
efficient. 

This work does not deal with ocean waves or wave and current combina- 
tions. Therefore, the effects of diffraction, free surface, or varying 
fluid velocity with depth along a vertical pile are not considered. 
Suffice it to note that the variation of the direction of the velocity 
vector at a given point with time and at a given time with depth tends 
to reduce the spanwise coherence along a vertical pile resulting in a 
reduced in-line as well as transverse force. Consequently, the trans- 
verse force coefficients presented herein represent their maximum 
Possible value since they have resulted from a highly correlated wake 


along the cylinder in a perfectly two-dimensional harmonic flow. 
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i ERIM ChSOR FORCESCOERFICIENTS 


In view of the fact that separated flow about bluff bodies is not 
yet amenable to theoretical analysis, it is necessary to obtain the force 
components experimentally through the use of appropriate coefficients. 
One can then demonstrate the dependence of these coefficients on the 
Parameters governing the flow and give some sense of universality to the 
data as well as the conclusions drawn from them. 

As pointed out earlier, the transverse force stems from the alter- 
nate shedding of vortices. During a given cycle of harmonic oscillation 
the velocity varies sinusoidally and hence the strength as well as the 
frequency of vortices. Thus, the time dependent transverse force con- 
tains first, second, and higher order harmonics. A very detailed study 
of the transverse force will require the evaluation of all the harmonics 
as well as their frequencies. Evidently this would be an extremely time-~ 
consuming process and would not necessarily add to the designer's ability 
to predict the forces acting on pilings to a greater degree of accuracy. 
In the present work only the maximum transverse force in a given cycle 
has been analyzed. A cursory examination of the force traces has shown 
that the second and higher order harmonics are considerably smaller. 
These second and higher order harmonics may be of some importance in the 
analysis of hydroelastic oscillations. 


The transverse force coefficient has been defined by: 


cLuay = Maximum peak of the transverse force in a cycle (1) 
(.5) (0) (D)(L) (U7) 
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The frequencies of oscillation have been expressed either in terms of 


the frequency of oscillation of the harmonic flow as 


— frequency of the first harmonic of the lift force in a cycle 
if frequency of the harmonic motion 


(2) 
or in terms of the Strouhal number defined by 
FD ; areal ; f 


ue nS (3) 


in which 9 represents the density of fluid; D and L, the diameter and 
length of the cylinder; Ua the amplitude of the oscillating velocity; 
and uve the frequency of the transverse force. 

As shown previously [9], CLMAX, St, and i depend on the Keulegan- 
Carpenter number defined by K = UT/Ds the Reynolds number defined by 
Re = UD/vs and the relative roughness defined by k/D. In the forego- 
ing, T represents the period of oscillation; v, the kinematic viscosity 
of the fluid; and k, the mean height of carefully-sieved sand grains. 


Thus one may write [9]: 


CLMAX 

UI 
ie = UF Sain a Re DO) (4) 
Sie 


It has also been shown previously [9] that either Re or UT/D in equa- 
tion (4) may be replaced by Re/K = Epo This parameter is called the 


frequency parameter and denoted by 8, so that 
ee 
B= D/vT E)) 


Evidently, for a series of experiments conducted with a cylinder of a 


given diameter in water (of uniform and constant temperature and density) 
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undergoing harmonic oscillations with a constant period (T), 8 is held 
constant. Then the variation of a force coefficient with K may be 
plotted for constant values of 8. Subsequently, one can easily recover 


the Reynolds number from 
Re = BK (6) 


and connect the points, on each 8 = constant curve, representing a given 
Reynolds number for a family of suitably selected values of the Reynolds 
number. Such a procedure eliminates the difficulty of trying to draw 
contours of constant K, or constant Re, or constant CLIMAX versus K or Re 
or K versus Re. Suffice it to note that the data reported herein shal] 


be analyzed according to the relationship 


kK 


[a coefficient] 


and the Reynolds number will be used in the manner cited above. The 
power of this new plotting procedure (new, as far as the wave force 


analysis is concerned) will become apparent later. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A. BRIEF DESCRIPTION OF THE EXPERIMENTAL EQUIPMENT 

The capabilities of a small U-shaped water tunnel have been demon- 
strated by Sarpkaya and Tuter [5] for generating oscillatory flow at 
relatively low Reynolds numbers. When it became desirable to investi- 
gate the in-line and transverse forces on cylinders at higher Reynolds 
numbers, a larger U-tunnel was constructed. 

The design and operational characteristics of this larger U-tunnel 
are described in great detail by Collins [9] and will not be repeated 
here in their entirety. A photograph of the tunnel used throughout this 
investigation is shown in Fig. 1. 

Of significance is the general shape of the tunnel and its effect 
upon the velocity of oscillation. The cross sectional area of the two 
legs (6x3 feet) is twice that of the test section (3x3 feet) so that the 
virtual amplitude” or the velocity of oscillation is exactly twice that 
of the free surface, Additionally, the length of the horizontal test 
Section is greater than twice the virtual amplitude so that fully de- 
veloped flow is achieved in the test section (see Fig. 2). 

Compressed air and four butterfly valves (operated by a pneumatic 
three-way control valve) were used to set the fluid into harmonic motion. 
It was considered prudent to allow the harmonic oscillations to damp out 
naturally during a test run rather than to employ the method of forced 
oscillations. The butterfly valves were operated only once, at the 


*Virtual amplitude is defined as the amplitude of oscillation which 
the cylinder experiences. 
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Photograph of the U-tunnel. 


Figure lI. 
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beginning of each run, causing a natural period of oscillation of 5.507 
seconds. The fluid oscillated smoothly, and the force and acceleration 
transducers required no filters whatsoever. 

Experiments were carried out at various water temperatures by heat- 
ing the fluid to the desired temperature in order to vary the Reynolds 


number for a given relative amplitude and relative surface roughness. 


B. CIRCULAR CYLINDER MODELS 

Circular cylinders with diameters ranging in size from 7 inches to 
2 inches were used in this investigation. The length of each cylinder 
was such that it allowed for a gap of 1/32 inch between the tunnel wall 
and each end of the cylinder. The cylinders were mounted and held in 
the test section by inserting the force transducers into each end of 
the cylinder. The force transducers used in this operation were identi- 
cal to those used in the previous study [9]. 

This investigation dealt solely with sand-roughened cylinders. The 
Sand was segregated by use of sieves of appropriate mesh sizes, so that 
the proper relative roughness, k/D, was achieved for each cylinder 
tested. This process resulted in uniform roughness which surpassed al] 
expectations. Photographs taken with a scanning electron microscope 
shox the uniformity of size and packing oF the sand grains (see Fig. 3). 


j 


C. FORCE MEASUREMENTS 

The instantaneous in-line and transverse forces were measured by two 
identical transducers. A special housing was fabricated for each gage 
So that it could be mounted on the tunnel window, inserted into the test 


cylinder, and rotated to measure either the in-line or the. transverse 


force, 
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Figure 3. 
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The gages had a nominal capacity of 500 Ibf with a 200 percent over- 
load capacity. Application of a 500 ibf load caused a deflection of .01 
inch. The largest actual load that the gages were subjected to during 
this investigation was less than 200 lbf, causing a deflection of .008 
inch. The gages were calibrated by suspending a load in the middle of 
the cylinder after setting each gage to sense only the transverse force. 
Repeated calibrations have shown the gages to be perfectly linear for 
the loads encountered in the investigation. 

Normally, one gage was used to measure the in-line force and the 
other to measure the transverse force, although at times both gages were 
used to measure only one. These force measurements were simultaneously 
recorded along with the instantaneous acceleration on two-channel Honey- 
well recorders operating at a speed of one inch per second. The natural 
period of 5.507 seconds allowed for 55.07. divisions per cycle. The ampli- 
tude of the transverse force and the flow characteristics such as K and 
Re were determined from these traces. Sample transverse force and ac- 


celeration traces are shown in Fig. 4. 


D. ACCELERATION MEASUREMENTS 

The method used during this investigation consisted of measuring the 
differential pressure between two pressure taps located as shown in Fig. 
5. These taps were placed symmetrically on the two vertical legs of the 
tunnel at an elevation of 50 inches below the mean water level so that 
H = 38 inches. By use of Bernoulli's equation, applied between each pres- 
sure tap and the instantaneous water level, it can be shown that the vir- 


tual amplitude (twice the amplitude of the free surface amplitude) is 


given by: 
A= 2A) = ——S (8) 
T- —(5)°H 
gq ul 
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Acceleration and transverse force traces. 


Figure 4, 
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Figure 5. Position of pressure taps. 


in which T and g are constant and H is kept constant. This method 
yields the virtual amplitude, or the maximum velocity in each cycle. 
The signal was entirely free of noise or small surface effects. The 
transducer was calibrated and its linearity checked before each series 
of experiments. 

Because this measurement is so critical, the accuracy of the method 
just described has been checked many times over. Each time the alter- 
nate method verified the accuracy of the procedure used. Iwo of these 
methods used for verification are described in [9]. A third check was 
provided by using a magnetic velocimeter which sensed changes in the 
Magnetic field caused by changes in the fluid velocity. As anticipated, 
the probe again verified the accuracy of the differential pressure 
method to within two percent accuracy, claimed by the manufacturer of 
the probe. 

These comparisons as well as the perfectly sinusoidal and noise free 
character of all pressure and force traces speak for the suitability of 


the unique test facility used in this study. 
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E. DATA REDUCTION 

Experiments were repeated at least twice for each cylinder, and nor- 
mally both of these runs were evaluated. Reynolds number, Keulegan- 
Carpenter number, and lift coefficient were calculated for between 20 
to 25 cycles for each run evaluated. 

It should be emphasized again that only the amplitude and frequency 


of the first harmonic of the transverse forces were calculated. 
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IV. DISCUSSION OF RESULTS 


The data obtained in the present investigation for various values 
of B and k/D are presented in Figs. 6 through 22 in terms of C, (max) 
versus K. A close examination of the data for a given k/D ratio and 8 
Shows that there is considerable scatter in the data due primarily to 
the random nature of the transverse force. Secondly, the variation of 
the transverse-force coefficient from one 8 to another, for a given 
cylinder and k/D ratio, is no more than the scatter in the data for a 
given 8. In other words, C, (max) does not depend, within the limits of 
accuracy of the data, on the Reynolds number particularly for relative 
rougnhnesses of k/D = 1/50, 1/100, and 1/200. The data obtained by 
Collins [9] also show that the transverse-force coefficient is indepen- 
dent even for a relative roughness as small as 1/400. 

The Reynolds-number independence of C, (max) for rough-walled cylin- 
ders, in the range of relative roughness from 1/400 to 1/50, is better 
demonstrated through the use of composite plots shown in Figs. 23 through 
25. Each figure contains al] data obtained with all cylinders for a 
given k/D. It is clearly apparent that the transverse-force coefficient 
1s independent of 8 and hence of the Reynolds number for the relative 
roughnesses indicated in these three figures. Figure 26 1s a similar 
plot for k/D = 1/800. It is apparent, within the limits of the scatter 
of the data, that there is some Reynolds number dependence, particularly 
for large values of 8. Additional data are needed at relative rough- 
nesses of 1/600, 1/1000, etc., in order to establish the variation of 


the transverse-force coefficient with Reynolds number as k/D approaches 
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zero. No attempt has been made in the present investigation to estab- 
lish the said relationship primarily because it is extremely difficult 
to find and apply very fine sand particles with sufficient accuracy on 
the cylinders used. From a practical viewpoint, relative roughnesses 
as small as 1/800 are not of special interest since the marine growth 
on the structural elements of offshore structures often give rise to 
considerably larger relative roughnesses (e.g. k/D = 100). 

One of the most remarkable effects of roughness on the transverse 
force coefficient is the increase of the said coefficient to values com- 
monly encountered at relatively low Reynolds numbers with smooth cylin- 
ders. This point is clearly demonstrated by plotting the mean tranverse 
force coefficient for smooth cylinders at relatively low 8 values on the 
lift-force traces shown in Fig. 25 for rough cylinders with k/D = 1/200. 
It should be noted in passing that any one of the three figures, namely 
Figs. 23-25, could have been used for this purpose since C, (max) does not 
vary significantly from one k/D to another for a given Keulegan-Carpenter 
number. Figure 25 shows that the transverse force coefficient for smooth 
cylinders at relatively low Reynolds numbers forms the upper limit of the 
Same coefficient for rough cylinders. This is an extremely important 
experimental finding particularly for model testing. It is evident from 
the foregoing that tests carried out at very low Reynolds numbers with 
smooth cylinders should yield information about the transverse force act- 
ing on rough-walled cylinders at very high Reynolds numbers. The exist- 
ence of such a relationship has been surmised but has never been clearly 
demonstrated as in the present investigation. 

As noted earlier, the alternating nature of the transverse force is 


as important as its magnitude. It is for this reason that the frequency 
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of the lift force has been determined for rough cylinders. A careful 
examination of the results have shown that IF gs or the Strouhal number 
remains essentially constant at a value of about 0.22. This result con- 
firms that arrived at independently by Collins [9] and will not be 
elaborated upon further. 

In considering the relevance of the coefficients presented herein 
to wave induced loads on offshore structures, it is of course important 
to take into account the differences between uniform two-dimensional 
harmonic motion and the wave motion where the velocity vector both rotates 
with time at a point and decays in magnitude with depth. The spanwise 
variations of the flow in general lead to reduced spanwise coherence. 
[t is safe to assume that both the three-dimensionality of the flow and 
the reduction of the correlation length along the cylinder, in an ocean 
environment, tend to increase the base pressure and thus give rise to 
force coefficients which are smaller than those obtained with purely two- 
dimensional flows. The transverse force coefficients presented herein 
obviously represent their maximum possible values since they have re- 
sulted from a uniform, two-dimensional flow where the instantaneous wake 
of the cylinder has the highest possible degree of Spanwise correlation. 
Thus, the value of the results presented herein lies in the fact that the 
designer now knows the maximum possible values of the transverse force if 
not the values which might be more appropriate to the conditions under 


which the structure must survive and function. 
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V. CONCLUSIONS 


The results presented herein warrant the following conclusions: 

(7) The transverse-force coefficient is nearly independent of 
the Reynolds number, particularly for relative roughnesses larger than 
about 1/800; 

(ii) The transverse force coefficient increases with increasing 
K in the range of K values from about 5 to 12 and then decreases with 
increasing K; 

(iii) The smooth cylinder data at relatively low values of the 
frequency parameter form more or less the upper limit of the rough 
cylinder data; 

(iv) The transverse force is a significant fraction of the total 
resistance at all Reynolds numbers and must be taken into considera- 
tion in the design of structures; 

(v) The Strouhal number remains nearly constant for all 
Reynolds numbers at about 0.22; 

(vi) The results presented herein and the conclusions arrived at 
are applicable only to cylinders in harmonic flow with zero mean veloc- 
ity within the range of Re, K, and k/D values encountered in the tests; and 

(vii) The force coefficients for wavy flows may differ somewhat 
from those presented herein partly due to the reduced coherence along 
the cylinder and partly due to the nonlinear interaction between the 


Current and the waves. 
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Figure 6. C,(max) versus K for k/D = 1/50, 3 inch. 
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Figure 7. C, (max) versus K for k/D = 1/50, 6 inch. 
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Figure 9. C, (max) versus K for k/D = 1/100, 2 1/2 inch. 
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C. (max) versus K for k/D = 1/100, 3 inch. 
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Figure 12. C, (max) versus K for k/D = 1/100, 5 inch. 
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Figure 13. C, (max) versus K for k/D = 1/100, 6 inch. 
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Figure 14. C, (max) versus K for k/D = 1/100, 7 inch. 
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Figure 15. C, (max) versus K for k/D = 1/200, 2 inch. 
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Figure 16. C, (max) versus K for k/D = 1/200, 2 1/2 inch. 
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Figure 17. C, (max) versus K for k/D = 1/200, 3 inch. 
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Figure 18. C, (max) versus K for k/D = 1/800, 3 inch. 
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Figure 19, C, (max) versus K for k/D = 1/800, 4 inch. 
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Figure 21. (C, (max) versus K for k/D = 1/800, 6 inch. 
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Figure 22. C, (max) versus K for k/D = 1/800, 7 inch. 
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